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The transverse Ising Model (TIM) in one dimension is the simplest model which exhibits a quan- 
tum phase transition (QPT). Quantities related to quantum information theoretic measures like 
entanglement, quantum discord (QD) and fidelity are known to provide signatures of QPTs. The 
issue is less well explored when the quantum system is subjected to decoherence due to its interac- 
tion, represented by a quantum channel, with an environment. In this paper we study the dynamics 
of the mutual information / (pab), the classical correlations C (pab) and the quantum correlations 
Q (pab), as measured by the QD, in a two-qubit state the density matrix of which is the reduced 
density matrix obtained from the ground state of the TIM in Id. The time evolution brought about 
by system-environment interactions is assumed to be Markovian in nature and the quantum channels 
considered are amplitude damping, bit-flip, phase-flip and bit-phase-flip. Each quantum channel is 
shown to be distinguished by a specific type of dynamics. In the case of the phase-flip channel, there 
is a finite time interval in which the quantum correlations are larger in magnitude than the classical 
correlations. For this channel as well as the bit-phase-flip channel, appropriate quantities associated 
with the dynamics of the correlations can be derived which signal the occurrence of a QPT. 
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I. INTRODUCTION 
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The correlations which exist between the different con- 

pH gtituents of an interacting quantum system have two dis- 

Qtinct components: classical and quantum. The most well- 

, I , jtnown example of quantum correlations is that of entan- 

r-| glement which serves as a fundamental resource in several 

d c(uantum information processing tasks [ll-Q- In the case 

? of bipartite quantum systems, the quantum discord (QD) 

I ^^j ias been proposed to quantify quantum correlations more 

general than those captured by entanglement [j-[y|. In fact, 

CSl there are separable mixed states which by definition are un- 

^ pntangled but have non-zero QD. The utility of such states 

O^ in certain computational tasks has recently been dcmon- 

^1^ strated both theoretically [3] and experimentally Q . QD 

-^ thus has the potential to serve as an important resource in 

• fcertain types of quantum information processing tasks. 

^__ . The quantum mutual information / (pab) measures the 

^•^ total correlations, with classical as well as quantum com- 

I ponents, in a bipartite quantum system and is given by 



I (pab) = S (pa) + S (pb) - S (pab) 



(1) 



p^ where pab is the density matrix of the full system and 
3 pA (pb) the reduced density matrix of subsystem A (B). 
■ - - Also, S{p) represents the von Neumann entropy with 
S{p) = -Tr{plog2p}. The QD , QipAB), is defined to 
be the difference between / (pab) and the classical corre- 
lations C (pab), i-e.. 
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C (pab) 



(2) 



The computation of classical correlations, C {pab)i is car- 
ried out in the following manner [3, M, llfl- I^ classical 
information theory, the mutual information / {A, B) = 
H{A) + H{B) - H{A,B) quantifies the total correla- 
tion between two random variables A and B. H{A) = 
-Ea-Palog2Pa, -^(-B) = -J2bPb^og2Pb and H{A,B) = 
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— ^^ ^ pa,b log2 Pa,b arc the Shannon entropies for the vari- 
able A, the variable B and the joint system AB respec- 
tively. The joint probability of the variables A and B 
having the values a and b respectively is represented by 
Pa,b and Pa = J2bPa,b^ Pb = J2aPa,b- The classical mu- 
tual information has an equivalent expression J{A, B) = 
H{A) — H{A\B) via the Bayes'rule. The conditional en- 
tropy H{A\B) is a measure of our ignorance about the 
state of A when that of B is known. In the case of a quan- 
tum system, the von Neumann entropy replaces the Shan- 
non entropy and the quantum generalization of the classi- 
cal mutual information I{A, B) is straightforward yielding 
the expression in equation ([1]). The quantum version of 
J {A, B) is not, however, equivalent to / (pab)- This is be- 
cause the magnitude of the quantum conditional entropy 
depends on the type of measurement performed on sub- 
system B to gain knowledge of its state so that different 
measurement choices yield different results. We consider 
von Neumann-type measurements on B defined in terms 
of a complete set of orthogonal projectors, 11^, correspond- 
ing to the set of possible outcomes i. The state of the 
system, once the measurement is made, is given by 



P^ 



{l<g,Uf)pAB{l<E>nf)/p, 



with 



p, = Tr((i®nf)p^B(/®nf)) 



(3) 



(4) 



Here / denotes the identity operator for the subsystem A 
and Pi gives the probability of obtaining the outcome i. 
The quantum analogue of the conditional entropy is 



S{pAB\{Uf})=J2p,S{p,) 



(5) 



The quantum extension of the classical mutual information 
is given by 



j{pAB\{Tlf})^SipA)^S{pAB\{Tlf}) 



(6) 



When projective measurements are made on the subsystem 
B, the non-classical correlations between the subsystems 
are removed. Since the value of J {pab\ {nf }) is depen- 



dent on the choice of {11^}, J should be maximized over ah 
{li-i} to ensure that it contains the whole of the classical 
correlations. Thus the quantity 



C(p^s)=inax(j(p^B|{nf})) 
{nf} 



(7) 



provides a quantitative measure of the total classical cor- 
relations @. 

Though the concept of the QD is firmly established, its 
computation is restricted to two-qubit states and that too 
when pAB has special forms [9l-[ll|. For two-qubit X-states, 
the density matrix in the basis {|11), |10), |01), |00)} has the 
general structure 



Px 



pu pi4 

P22 P23 

P32 P33 

P41 P44 



(8) 



with pi2 = P21 = Pl3 = P31 = P24 = P42 = P34 = P43 = 0. 

Analytic expressions for the QD can be derived only in 
some special cases. We restrict attention to the case 



PAB 



a / 
b z 
z b 

food. 



(9) 



where A, B correspond to the two individual qubits and z, 
/ are real numbers. The eigenvalues of pab are [l3| 

Ao = ^ I (1 + C3) + \IacI + {^^^^^'' 



with 



4 

A2 = T (1 - C3 + Cl +C2) 
A3 = - (1 - C3 - Ci - C2) 



ci = 22 + 2/ 
C2 = 2z- 2/ 
C3 = a + (i — 25 
C4 = a — d 



(10) 



(11) 



The mutual information / {pab) (equation ([I}) can be writ- 
ten as d, [3 



/ {pab) = S (pa) + S (pb) + V A„ log2 A„ (12) 



a=0 



where 



S{pa) ^ S{pb) 



1 -f C4 1 -I- C4 

1 — C4 1 — C4 

— ;;— log2^:;— 



(13) 



With the expressions for I {pab) and C {pab) given in 
equations ([T^ . ([T^ and (O respectively, the QD, Q {pab), 
(equation 1^) can in principle be computed. The difficulty 
lies in the maximization procedure to be carried out in or- 
der to compute C {pab)- When pab is of the form given 



in ([H]), the maximization can be done analytically |12| re- 
sulting in the following expression for the QD; 



Q{pab) = min{(3i,(32} 



(14) 



where 



Qi = S {pb) - S {pab) - alog2 — — -&log2 — -r 

a + b a + 
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and 



Q2 = S {pb) - S{pab) 
-A_ log2 A_ 



A+log 



(15) 



(16) 



with A± = i (1 ± r) and T^ = (a - df + 4 (|z| + \f\f 

Quantum systems, in general, are open systems because 
of the inevitable interaction of a system with its environ- 
ment. This results in decoherence, i.e., a gradual loss from 
a coherent superposition to a statistical mixture with an 
accompanying decay of the quantum correlations in com- 
posite systems. The dynamics of entanglement and QD 
under system-environment interactions have been investi- 
gated in a number of recent studies J13l - [la |. One feature 
which emerges out of such studies is that the QD is more 
robust than entanglement in the case of Markovian (mcm- 
oryless) dynamics. The dynamics may bring about the 
complete disappearance of entanglement at a finite time 
termed the 'entanglement sudden death' |13l. Il4| . The QD, 
however, is found to decay in time but vanishes only asymp- 
totically fla. [l7t - ll9l | . Also, under Markovian time evolution 
and for a class of states, the decay rates of the classical and 
quantum correlations exhibit sudden changes in behaviour 
[16l . |l7[. Three general types of dynamics under the ef- 
fect of decoherence have been observed [ia|: (i) C {pab) 
is constant and Q {pab) decays monotonically as a func- 
tion of time, (ii) C {pab) decays monotonically over time 
till a parametrized time p^c and remains constant there- 
after. Q {pab) has an abrupt change in the decay rate 
at psc and has magnitude greater than that of C {pab) in 
a parametrized time interval and (iii) both C {pab) and 
Q {pab) decay monotonically. Mazzola et al. [13| have 
demonstrated that under Markovian dynamics (qubits in- 
teracting with non-dissipative independent reservoirs) and 
for a class of initial states the QD remains constant in a 
finite time interval < t < t. In this time interval, the 
classical correlations, C{pab), decay monotonically. Be- 
yond t = i, C {pab) becomes constant while the QD de- 
creases monotonically with time. The sudden change in 
the decay rate of correlations and their constancy in cer- 
tain time intervals have been demonstrated in actual ex- 
periments [23,[2l|. 

In this paper, we focus on a two-qubit system each qubit 
of which interacts with an independent reservoir. The den- 
sity matrix of the two-qubit system is described by the re- 
duced density matrix derived from the ground state density 
matrix of the transverse Ising model (TIM) in one dimen- 
sion (Id). We investigate the dynamics of the QD as well as 
the classical correlations under Markovian time evolution 
and identify some new features close to the quantum criti- 
cal point of the model. In Sec. II, the calculational scheme 
for studying the dynamics of the classical and quantum 
correlations is introduced. We further describe the quan- 



turn channels representing the system-environment inter- 
actions for which the computations are carried out. Sec. 
Ill presents the major resuhs obtained and a discussion 
thereof. Sec. IV contains some concluding remarks. 



II. DYNAMICS OF CLASSICAL AND QUANTUM 
CORRELATIONS 



We consider the TIM Hamiltonian in Id described by 
the Hamiltonian 



0.18 



H 






L 



(17) 



where a^ and ct^ are the Pauli matrices defined at the site 
J of the chain and L is the total number of sites in the 
chain. We further assume periodic boundary conditions. 
The Hamiltonian in (|17p can be exactly diagonalized in 
the thermodynamic limit L — >■ oo [22, [23] . When the pa- 
rameter A = 0, all the spins are oriented in the positive 
z direction in the ground state whereas for A = oo, the 
ground state is doubly degenerate with all the spin point- 
ing in either the positive or the negative direction. As one 
goes from one limit to the other, a quantum phase transi- 
tion (QPT) occurs at the critical point Ac = 1 separating 
two different phases, the paramagnetic phase with the mag- 
netization (cr^) zero and the ordered ferromagnetic phase 
characterized by a non-zero magnetization. The QPT is 
signaled by the divergence of the correlation length at the 
critical point. Since the ground state wave function un- 
dergoes a qualitative change at the critical point, it is rea- 
sonable to expect that the quantum correlations present 
in the ground state would provide signatures of the oc- 
currence of a QPT. Such signatures in fact do exist for 
different measures of quantum correlations, namely, entan- 
glement [H, 0, Hi^il and QD [I^, HI [Ig- For the TIM, 
the two-site reduced densit y rn atrix p^ has the form given 
in equation ® [H, [H, [H [13 with 
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The magnetization {a") of the TIM is given by |22l.[25l 
1 r .,(l-l-Acos 



(18) 



(a^) 



W0 



where 



u;0 = Y (Asinc/)) +(1 + A cost 



(19) 



(20) 



is the energy spectrum. The spin-spin correlation functions 
are obtained from the determinant of Toeplitz matrices J25l . 




Figure 1. (Color online) Amplitude damping channel: decay 
of mutual information / [pab) (solid line), classical correlations 
C (pab) (dashed line) and quantum correlations (QD) Q (pab) 
(dotted line) as a function of parametrized time p = 1 - 
and A = 0.5 
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(21) 
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We next consider the interaction of the chain of qubits, 
each qubit representing an Ising spin, with an environment. 
We choose the initial state of the whole system at time 
i = to be of the product form, i.e.. 



p{0) = Ps{0) (S PeiO) 



(23) 



where the density matrices ps and pe correspond to the 
system and environment respectively. We assume that the 
environment is represented by L independent reservoirs 
each of which interacts locally with a qubit constituting 
the Ising chain. The two-qubit reduced density matrix ob- 
tained from equation (|23p can be written as 

(24) 



Pr(0) = prs{0) ^ PreiO) 



where prs and pre represent the two-qubit reduced density 
matrix of the transverse Ising chain and the corresponding 
reduced density matrix of the two-reservoir environment 
respectively. The two-qubit reduced density matrix prs is 
obtained by taking partial trace on ps over the states of all 
the qubits other than the two chosen qubits. Similarly, p^e 
is obtained from pe by taking a partial trace over the states 



of all the reservoirs other than the two local reservoirs of 
the selected qubits. The quantum channel describing the 
interaction between a qubit and its environment can be 
of various types: amplitude damping, phase damping, bit- 
flip, phase-flip, bit-phase-flip etc. [1^[33- Our objective is 
to investigate the dynamics of the two-qubit classical and 
quantum correlations (in the form of the QD) under the 
influence of various quantum channels. 

The time evolution of the closed quantum system, com- 
prised of both the system and the environment, is given 

by 



Pseit)=Uit)p,,{0)U^{t) 



(25) 



where U{t) is the unitary evolution operator generated by 
the total Hamihonian H of the system (U = e"*^*/'^). H 
is given hy H = Hg + H^ + Hse where Hg and He rep- 
resent the bare system and environment Hamiltonians re- 
spectively and Hse the Hamiltonian describing the interac- 
tions between the system and the environment. The time 
evolution of the system s subject to the influence of the 
environment e is obtained by carrying out a partial trace 
on Pse{t) (equation (|25p ) over the environment states, i.e.. 



pg{t) = Tre [U{t)pge{m\i)\ 



(26) 



Let \ek) be an orthogonal basis spanning the finite- 
dimensional state space of the environment. With the ini- 
tial state of the whole system given by equation 



pS) = Jl^^fcl^ [/'-('^) ® /'<=(o)] ^^l^fc) 



(27) 



Let /9e(0) = |eo)(eo| be the initial state of the environment. 
Then 



p,(t) = 5^i?,p.(0)4 



(28) 



where Ek = (efe|L'^|eo) is the Kraus operator which acts on 
the state space of the system only [ij, [SO]- Let {0i}, i = 
1, 2, ..., d, define the basis in the state space of the system 
s. There are then at most cP independent Kraus operators 
Ek, fc = 0, ..., d^ - 1 [30,[m. The unitary evolution of s + e 
is given by the map: 



\^i)\eo) -> Eo\<f>i)\eo) 
|</'2)|eo)->So|02)|eo) 



Ed^-i\4>i)\ed2-i) 
Ed^-i\(l)2)\ed^-i) 



|0d)|eo> ^ Eo\(f>d)\eo) + - + Ed2_i\(j>d)\ed2^i) (29) 
In compact notation, the map is given by 

C/|0,)|eo)=^Sfe|0.)|e;c), i = l,2,...,d (30) 

k 

In the case of N system parts with each part interacting 
with a local independent environment, equation (|28p be- 
comes 
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(31) 



where Ej°' is the kath Kraus operator with the environ- 
ment acting on system part a. The specific form for Ps{t) 




Figure 2. (Color online) Bit-flip channel: decay of mutual 
information I (pab) (solid line), classical correlations C (pab) 
(dashed line) and quantum correlations (QD) Q {pab) (dotted 
line) as a function of parametrized timep = 1 — e"''* and A = 0.5 



arises as the total evolution operator can be written as 
U{t) = Ui{t) (» U2{t) (» ... ® UN{t). Following the gen- 
eral formalism of the Kraus operator representation, an 
initial state, Prs(O), of the two-qubit reduced density ma- 
trix evolves as [19, ^ 



,{t)=J2E^,r.PrsiO)El, 



(32) 



flM 



where the Kraus operators Efj^,^ = E^ ® Ei, satisfy the 
completeness relation ^ ^ E^^^EJ^ u = I for all t. We now 
briefly describe the various quantum channels considered 
in the paper and write down the corresponding Kraus op- 
erators. A fuller description can be obtained from Refs. 

[il,|3^. 

(i) Amplitude Damping Channel. 

The channel describes the dissipative interaction be- 
tween a system and its environment resulting in an ex- 
change of energy between s and e so that s is ultimately 
in thermal equilibrium with e. The s + e time evolution is 
given by the unitary transformation 



|0),|0)e^|0),|0). 



Vg|l).,|0)e + VP|0)s|i) 



(33) 



(34) 



where |0)s and |l)s are the ground and excited qubit states 
and |0)e, |l)e denote states of the environment with no 
excitation (vacuum state) and one excitation respectively. 
Equation (j33p stipulates that there is no dynamic evolu- 
tion if the system and the environment are in their ground 
states. Equation ([M|) states that if the system qubit is in 
the excited state, the probability to remain in the same 
state is q and the probability of decaying to the ground 
state is p [p + q = 1). The decay of the qubit state is 
accompanied by a transition of the environment to a state 
with one excitation. The qubit states may be two atomic 
states with the excited state decaying to the ground state 
by emitting a photon. The environment on acquiring the 
photon is no longer in the vacuum state. With a knowl- 
edge of the map equations (equations (p3)) and (p4| '). the 
Kraus operators for the amplitude damping channel can be 



0.18 





Figure 3. (Color online) Phase-flip channel: decay of mutual 
information I {pab) (solid line), classical correlations C (pab) 
(dashed line) and quantum correlations (QD) Q (pab) (dotted 



line) as a function of parametrized time p = 1 
0.5. Also, p„i = 0.0932, p^c = 0.1347 and p^ra 



-7* 



and A 



: 0.1649 



Figure 4. (Color online) Bit-phase-flip channel: decay of mutual 
information I (pab) (solid line), classical correlations C {pab) 
(dashed line) and quantum correlations (QD) Q [pab) (dotted 
line) as a function of parametrized time p = 1 
0.5. Also, psc = 0.0666 
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written as 



Bit-phase-flip 



Er,= 
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(35) 



where q ~ 1 — p. The Kraus operators for the two distinct 
environments (one for each qubit) have identical forms. In 
the case of Markovian time evolution, p is given by p — 
1 — e"'''* with 7 denoting the decay rate. 

(ii) Phase Damping (dephasing) Channel. The channel de- 
scribes the loss of quantum coherence without loss of en- 
ergy. The Kraus operators are: 



£^0 = 



1 

V9 



Ei = 





VP 



(36) 



with q= 1 — p and p = 1 — e '''*. 

(iii) Bit-flip, phase-flip and bit-phase-flip channels. The 
channels destroy the information contained in the phase 
relations without involving an exchange of energy. The 
Kraus operators are 



£n = 



1 
1 



Ei^ 



(37) 



where i = x for the bit-flip, i — y for the bit-phase-flip 
and i = z for the phase-flip channel with q' = 1 — p/2 and 
p = \ — e"'*'*. The expanded forms of the Kraus operators 
are: 

Bit-flip 



Phase- flip 
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(39) 



v/T^p72 

/r^W2 

-i\fpj2' 




(40) 



As shown in Ref. [30|, the phase damping quantum oper- 
ation is identical to that of the phase-flip channel so that 
we will consider only one of these, the phase-flip channel, 
in the following. 

For a specific quantum channel, it is now straightforward 
to calculate the dynamics of the classical and quantum cor- 
relations. Equation (^^ describes the time evolution of the 
reduced density matrix of the TIM subjected to the influ- 
ence of an environment via a quantum channel. The initial 
state prs (0) has the form given in equation ^^ the elements 
of which are known via the equations (|18p - (P^ . The time- 
evolved state Prs{t) has again the form given in equation 
([9]) with the time dependence occurring in only the off- 
diagonal elements. With a knowledge of the elements, the 
mutual information I (pab), the QD Q{pab) and the clas- 
sical correlations C{pab) can be computed at any time t 
with the help of the formulae in equations ([TU| - p^ and 
equation ^. The results of our calculations for the vari- 
ous quantum channels are described in the next Section. 



III. RESULTS 

In the following, we make the substitution pab — Prsit). 
Amplitude Damping Channel. The dynamical evolution 
of the mutual information I (pab), the classical correla- 
tions C{pab) and the QD Q{pab) as a function of the 
parametrized time p {p = I — e^^*") is shown in FiglT] for 
A = 0.5. The solid, dashed and dotted lines represent 
the variations oi I[pab), C{pab) and Q{pab) respectively 
with p. All the correlations decay to zero in the asymp- 
totic limit of t — >■ oo, i.e., p — >■ 1. There is further no 
parametrized time interval or point when the quantum cor- 
relation becomes greater than the classical correlation. 





Figure 5. (Color online) Variations of pcri (solid line) and Pcr2 
(dashed line) with A in the case of the phase-flip channel; (inset) 
The first derivative of pcri and pcr2 w.r.t A diverges as the QCP 
Ac = 1 is approached 



Figure 6. (Color online) Variations of psc with A for the phase- 
flip (dashed line) and the bit-phase-flip (solid line) channels; 
(inset) The first derivative of psc w.r.t A exhibits a divergence 
as the QCP Ac = 1 is approached 
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Bit-flip Channel. Figj2] exhibits the dynamical evolution 
of I{pab) (solid line), C{pab) (dashed line) and Q{pab) 
(dotted line) as a function of the parametrized times and 
with A = 0.5. The quantum correlations disappear com- 
pletely in the asymptotic limit p ^ 1. In the same limit, 
I{pab) = C{pab) has a finite value. In the case of both the 
amplitude damping and bit-flip channels, the same features 
as observed for A = 0.5 are obtained for the other values 
of A. 

Phase-flip Channel. In Fig|31 we plot the variations of 
I (pab) (sohd line), C (pab) (dashed line) and Q (pab) 
(dotted line) versus the parametrized time p with A = 0.5. 
There is a sudden change in the decay rates of both C (pab) 
and Q (pab) a.t p = Psc. There are two points, p = Pcn and 
P — Pcr2 at which the plots of C (pab) and Q (pab) cross 
each other. The classical correlations remains constant be- 
yond the point p = Psc whereas the QD decays asymptoti- 
cally to zero. In the parametrized time interval Psc < P < 
1, the magnitude of C {pab) = -^(pa_b)L=i, the mutual 
information of the completely decohered state {p =1). In 
the interval Pcn < p < Pcr^ ; the quantum correlations are 
larger in magnitude than the classical correlations contra- 
dicting an earlier conjecture that C (pab) > Q (pab) in 
any quantum state [l6|. At the crossing points, Pcn and 
Pcr2^ one gets the equality C {pab) = Q {p/ "^ ~ -^° 



wcr2i 

Xu et al 



^AB 



20[ have recently investigated the dynamics of 
classical and quantum correlations under decoherence in 
an all-optical experimental setup. Fig. 4 of their paper pro- 
vides experimental verification of the dynamics displayed 
in FiglSl 

The sudden change in the decay rates of Q {pab ) and 
C {pab) at p = Psc is understood by noting that for p < 
Psc, Q = Q2 (equation dH])) and for p > psc, Q = Qi 
with C {pab) given by C {pab) = I {pab) - Q {pab)- At 
the crossing points Pcn and Pcr2, Q{pab) = C {pab) so 
that Pen and PcT-2 are the solutions of the equations Q2 = 
{I {pab) - Q2) and Qi = (/ {pab) - Qi) respectively. The 
constancy of C {pab) for values of p > Psc is explained by 
the fact that Q = Qi in this regime. From equations ([1]), 



C (pab) = S {pa) + a log2 ■ 

+ rfl0g2 



^fog2— TT 
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blog2 



d + b 



(41) 



As already pointed out, the time-evolved state Prs{t) has 
the form given in equation (jH]) with the diagonal elements 
a, b and d being independent of time. From (|lip and (jl3|) . 
S (pb) is thus independent of time. The other terms in 
equation (|4ip are also independent of time since they in- 
volve only the elements a, b and d. 

Bit-phase-flip Channel. Fig [3] shows the plots of I (pab) 
(solid hue), C {pab) (dashed line) and Q {pab) (dotted 
line) as a function of p. Again, there is a sudden change, as 
in the case of the phase flip channel, in the decay dynamics 
of C (pab) and Q (pab) at p = Psc but in this case the two 
plots do not cross each other but touch at a single point 

P = Psc- 

The dynamical features for the different quantum chan- 
nels have been reported earlier [la| for the class of states 
with a = d in the reduced density matrix of equation 1^ , 
i.e., C4 = in equation (jlip . In our study, the reduced den- 
sity matrix has the form shown in equation ^ with a ^ d, 
i.e., C4 7^ 0. The form corresponds to that of the two-qubit 
reduced density matrix obtained from the ground state of 
the TIM in Id. In this case each quantum channel is dis- 
tinguished by a specific type of dynamics. In Ref. [I3, the 
parameters ci , C2 and C3 are free and different types of dy- 
namics occur in different parameter regions corresponding 
to the same quantum channel. 

We now present some totally new results which have not 
been reported earlier. FiglS] shows plots oi Pcn (solid line) 
and Pcr2 (dashed line) versus A, the parameter appearing 
in the TIM Hamiltonian (equation p7|) ). in the case of the 
phase-flip channel. The inset of the Figure shows that the 
flrst derivative of Pcn and Pcr2 (both of which depend on 
A) w.r.t. the parameter A diverges as the QCP Ac = 1 is 
approached. The observation identifles a quantity which 
provides the signature of a QPT occurring in a system 
subjected to decoherence under Markovian time evolution. 
Fig[5] shows a variation of p^c with A for the phase-flip 
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Figure 7. (Color online) Variation of Apcr = Pcr2 — Pcri with 
A; (inset) The first derivative of Apcr w.r.t A diverges in the 
negative direction as the QCP Ac = 1 is approached 



(dashed line) and the bit-phase-flip (solid line) channels. 
Again, the inset shows that the first derivative of Psc w.r.t 
A diverges as the QCP Ac = 1 is approached. Figl7]exhibits 
the plot of Apcr = Pcr2 ~ Pcri versus A in the case of the 
phase-flip channel. The inset shows that ^ f^^''' diverges 
in the negative direction as the QCP Ac = 1 is approached. 
We remind ourselves that when p falls in the interval Apcr , 
the quantum correlations Q (pab) ai'c larger in magnitude 
than the classical correlations C{pab)- Figgis an out- 
conic of the results of Figl5]as Apcr = Pcr2 ~ Pcrx ■ In sum- 
mary, the first derivative of any one of the quantities Pcn , 
Pcr2 1 Psc and Apcr w.r.t A signals a quantum critical point 
transition. The quantities correspond to states in which 
the quantum correlations are either equal to or greater 
than the classical correlations. At p — pcn and p = Pcr2i 
Q{Pab) = C (pab) = I (pab) /2 which is characteristic 
of pure states with Q (pab) being equal to the entropy of 
entanglement [5|, |32| . The amplitude damping and bit- flip 
channels do not have these features. The appearance of 
singularities in the first derivatives of the quantities Pcn , 
Pcr2 1 Psc and Apcr with respect to the tuning parameter as 
the quantum critical point is approached can be explained 
by the fact that these quantities depend on two-spin cor- 
relation functions which exhibit a similar property close to 



the critical point. The non-trivial aspect arises from the 
identification of appropriate quantities associated with the 
dynamics of correlations which provide clear signatures of 
QPTs. 



IV. CONCLUDING REMARKS 

The TIM in Id is a prototypical example of a quantum 
system exhibiting a QPT. The many body ground state has 
both classical and quantum correlations. The QD provides 
a quantitative measure of the quantum correlations in a 
two-qubit state. In this paper, we consider a two-qubit 
state described by the reduced density matrix obtained 
from the ground state of the TIM in Id. The two-qubit 
state undergoes Markovian time evolution, described by 
the Kraus operator formalism, due to the local interactions 
of the qubits with independent environments. We consider 
the quantum channels, representing the interactions, to be 
of four types: amplitude damping, bit-flip, phase-flip and 
bit-phase- flip. The dynamics of the classical and quantum 
correlations exhibit distinctive features for each quantum 
channel. These features have been reported in an earlier 
study [1^ for a different class of initial states. In our study, 
we have not found evidence of another type of dynamical 
behaviour mentioned in [l6[ , namely, that the classical cor- 
relations, C (pab), are independent of time throughout the 
parametrized time interval p whereas the QD, Q (pab), de- 
creases monotonically and becomes zero in the asymptotic 
limit p — >■ I. The time evolution of the reduced density 
matrix of the TIM in Id further docs not exhibit the inter- 
esting dynamical behaviour described in [l7|, namely, the 
existence of intervals of parametrized time when C (pab) 
and Q (pab) are individually frozen. The most significant 
result of our study is the identification of quantities as- 
sociated with the dynamics of the classical and quantum 
correlations which diverge as the QCP of the TIM in Id, 
Ac = 1 , is approached thus providing a distinctive signature 
of a QPT from a different perspective. The generalization 
of the result to other model systems exhibiting QPTs will 
certainly be of significant interest. In the present study, we 
have restricted our attention to Markovian time evolution. 
The more general case of non-Markovian time evolution 
can be investigated only after the appropriate calculational 
scheme for a system of interacting qubits is developed [l^l • 
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